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ABSTRACT 


A previously developed code for calculating the mobility of 
charge carriers in narrow bandgap semiconductors does not predict 
the correct temperature dependence in all cases. It is thought 
that this is due to the way the electronic screening of the 
carriers is treated in the model. The objective of this research 
is to improve the handling of the screening by going beyond the 
current first Born approximation. Much of this work is directly 
related to the alloy semiconductor Hg. Cd Te wh i° h is important 
for infrared detectors and is a good candidate for microgravity 
crystal growth. The principal conclusion, so far, is that the 
major difficulty is probably the treatment of short range screen- 
ing at higher temperatures. 
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1 . 


Introduction 


H< ?-. x Cd Te is a mixed pseudo-binary semiconductor. Its 
energy gap fcafi Kfe "tuned" from 0.0 to 1.6 eV as x varies from x 
.17 to x=l at zero degrees (Dornhaus, Nimtz and B. Schlicht, 
1983, p. 148). This variation of band gap with x allows one to 
optimize for infrared detection (Long and Schmit, 1970, Kruse, 
1981) . Because of difficulties with producing compositionally 
uniform crystals due to gravity driven convective effects, 
H 9 n_ x Cd x Te is a good candidate for production in the reduced 
gravity environment of space (Lehoczky and Szofran, 1981, 1982). 


Semiconductors with small gaps tend to have small 
effective masses and hence high mobilities (Kruse, 1981) . The 
band structure of these materials is now understood from Kane's 
three level band model and k.p perturbation calculation (Kane, 
1957) which leads to non parabolic band structure. See appendix 
A. 5. 


Since these semiconductors are polar, scattering by the 
optical modes is important. The energy of the optical phonons is 
comparable to the energy of the electrons and hence the electrons 
are not scattered elastically. This means a relaxation time 
approximation is not valid (Howard and Sondheimer, 1953; Dingle, 
1956) and therefore variational methods are used to solve the 
Boltzmann equation (Ehrenreich, 1957) . The Boltzmann equation 
must contain all appropriate scattering terms and the narrow band 
gap means that thermally excited electrons will screen the 
carriers in certain scattering interactions (Ehrenreich, 1959 - 
two papers) . We attempt to improve the calculation of the 
scattering by improving the treatment of screening (see Appendix 
A. 1) . There are several contributions to the scattering of 
electrons in Hg, Cd Te and hence to the determination of 
electron mobility (Xeh<bczky, Szofran and Martin, 1980; Lehoczky, 
Summers, Szofran and Martin, 1982). These include longitudinal - 
optic phonon interactions (LO) , acoustic phonon scattering (ac) , 
ionized impurity scattering (ii) , electron hole scattering (eh), 
compositional disorder scattering (dis) and neutral defect (nd) 
scattering. Screening is involved in (LO) , (ii) , and (eh) . At 
temperatures above approximately liquid nitrogen temperatures, LO 
scattering is easily the most important. See Appendix A. 6 for a 
discussion of the basic interaction term in the Hamiltonian. The 
LO modes are important rather than the TO modes because it is 
only the LO modes that have strong electric fields which 
accompany their vibration 
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2 . 


In alloy semiconductors, such as Hg._ x Cd Te, there is a 
strong polar coupling of the carriers to the optic modes. In 
describing this interaction, one must properly account for the 
screening by use of a suitable dielectric function (Ziman, 1972; 
Mahan, 1981; Whitsett et al., 1981, see also appendix A.2) . The 
dielectric function depends on several factors (Lehoczky et a . 
1974 Nelson et al, 1978) among which is the Fermi Thomas wave 
vector or its reciprocal which is essentially the screening 
length. Previous calculations of electron mobility considering, 
among other factors, the scattering of carriers by optical Phonon 
modes have yielded electron mobilities which are too high 
(Lehoczky et al, 1982). We aim to improve these calculations by 
using the Friedel sum rule and the second Born approximation 
(Joachain, 1975) for the phase shifts to calculate a better 
screening length (Stern, 1967; Agarwal and Singh, 1983). 


A second objective is to fill in 
ground material by use of a suitable set 
reviews are by Zawadski, 1974; Chattopadhyay 
Chattopadhyay and Nag, 1974; Nag, 1980; Nag, 


the necessary back- 
of appendices (some 
and Queisser, 1981; 
1975; Rode, 1975) . 


It is essential in all our calculations to take into 
account the fact that the energy bands are non parabolic, 
although they may still be assumed to be spherically symmetric 
(Szymanska and Dietl, 1978; Pfeffer and Zawadzki, 1978) . The non 
parabolicity of the band is a major complicating feature. 
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3. 


The Calculation 


The background details and the definitions can be found 
in the appendices. Here we outline the calculation. The idea is 
to choose a form for the screened potential of a charged impurity 
and then to use the Friedel sum rule to place a constraint on 
this potential which involves the screening length. From 
Appendix A. 3, the Friedel sum rule is 




*/r 2L 




(0 


From Appendix A. 4 the phase shifts // can be determined, at least 
to second order Born approximation accuracy, by 


where 



kJU 

I - 


A* 
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with U(r) = 



rh*/*’) V(t), 


(4b) 


and V(r) is 

Appendix A.l) 

V(r) = 


the customary screened Coulomb potential (see 

- Za e r/L / r , (5) 


with L being the screening length. 
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Gu,t') = (ter 7 ). 


with f<, Y> being the smaller and the larger of f and f'and 
are the spherical Bessel functions. 

Since ~ta / k = through second order we can write Eq. 

(2) as 



1 — If] £ J r 


Since 


2Jj( (2 Ji -m) x ) = i. 


we can easily show 

S x (2^-t-O - - ZtK L . 


( 6 ) 


(7) 


If we assume (4r) «1 for r4 L where L is the screening 
length in the spherical Bessel functions of Eqs. (3) and (4) (see 
Appendix A. 7 for a discussion of the validity of our approxima- 
tions) then we find 
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Following Agarwal and Singh 
A. 7), we approximate Eq. (6) as 


(1983) (also see Appendix 
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The screening length 
where 


in the first 


Born approximation is 







l = */> *,* 


( 11 ) 


(see e.q. Szymanska and Dietl, 1978) This expression is good even 
for non parabolic bands (See Appendix A. 8 for different ways of 
expressing the density of states) . We find 

( L/u ,)" = / - , 02 ) 


where m* has been treated as constant in the factor multiplying 
in Eq. (2) . 

Eq. (12) is readily solved with the following results 

(l/ Lm ) = fiT7'‘-4, 


a I- 


which gives with A = _ j- 


L 


the following numbers. 


A 

VL (() 

-1 

2.41 

-.5 

1.62 

0 

1 

.5 

. 62 

1 

.41 
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4. 


Conclusions and Recommendations 


The results of the previous section were derived with 
the following approximations applied to the correction term to 
the first Born approximation [ /( t ~ 3 e /rf.) of Ec 3* ( 6 ) 3 

a. A constant effective mass can be used beyond 
the first Born approximation, 

b. -kL«^f where L is the screening length 

c. need only the^ = 0 term. 

We have discussed this approximation in A. 7. In Appendix A. 9 we 
note that improvements are not so easy and it is not clear our 
results are valid at higher temperatures. 

The major comment to make about this result is that it 
disagrees with that of Agarwal and Singh (19 83) who in my 
notation get 

(L/U.f = 


Both Eq. (12) (my result) and Eq. (13) assume lLdZl<^L 
It should be noted that my result differs qualitatively from 
Agarwal and Singh. They predict corrections to the first Born 
approximation increase L (for donors, Z>0) whereas I predict they 
decrease L. The literature seems to be divided as to which 
qualitative effect to expect (D. Chattopadhyay and H.J. Queisser, 
1981) . It should be mentioned that a decrease in L would mean 
more screening, less scattering and hence higher mobility - the 
opposite result from what we expected. 

We can also obtain the results of Agarwal and Singh if 
we use their expression for the second Born approximation but 
this expression does not appear to agree with standard sources 
(Joachain, 1975) . 

It remains to be seen whether better treatment of 
corrections to the first Born approximation would lower the 
mobilities and yield better agreement with experiment. This 
would be the first thing we would recommend investigating as our 
low energy approximation are suspect at higher temperatures. As 
shown in Appendix A. 9, going beyond M = 0 in the correction term 
quickly yielded equations that must be solved numerically for L. 


-/ 


( / — L ?-zJ 
/ ~f~ D Z: t 


(13) 
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Also if we drop the assumption that&L«l, then the integrals may 
still be doable, but they certainly are not convenient. 

Experimentally the calculations of the mobility are off 
by a factor of perhaps 2. If the corrections to the first Born 
approximation were enough to bring agreement with experiment it 
is likely that we would have cause to suspect the convergence of 
our procedure. 

Our recommendations would stress the positive aspects of 
the calculation which has already been done. The calculation is 
so complex that it is remarkable that qualitative and in some 
cases perhaps quantitative agreement with experiment can be 
achieved. It would appear that a review of the whole calculation 
is in order. It could be that the lack of agreement with experi- 
ment might not reflect just one error but perhaps several, none 
of them particularly serious by themselves. 

The error at high temperature should mean that close in 
scattering is overscreened by this model. Whatever the approach 
it seems this basic problem should be fixed. We also need to 
take into account that in general more than one type of carrier 
can contribute to electric current and to screening. There is 
also the troublesome problem of antiscreening which can occur in 
polar optic mode scattering when the electrons don't move fast 
enough to "follow the charge motion due to lattice waves 
(Ehrenreich, 1959) . 
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Appendix A. 1 
Screening 


We present here the basic ideas (Ziman, 1972) without 
putting in the details. A key idea we wish to mention is the 
difference between screening by free electrons and by band 
electrons which are described by Bloch Wave functions. 


For a free electron gas (with a uniform positive background) 
the applied potential can be written 


V(T\t) 




.t) 


d cl cO 


(A. 1 . 1 ) 


with r being the space coordinate, t the time, cf and u> are the 
usual wave vector and frequency.^ The Fourier components of the 
effective potential is then V(q,to) divided by the dielectric 
function. 


The Lindhard or Random Phase Approximation dielectric 
function is the Fourier component of the displacement vector 
divided by the Fourier component of the electric field and it 
comes out; 


e(v-). 1 - 


mre z 







W> t- ^ 

k * ‘j 


(A.1 .2) 


where f is the Fermi Function and Er> is the electron energy at 
wave vector ^ . * 

It is particularly ^interesting to examine this for the 
static case (CO = 0) when q is near zero. For then (assuming one 
band) 


E r - s - f ■ ^ * 

_ - / H 

4n?i- - ~ V 
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The static dielectric constant then comes out 

£ 


£(f,o) - I +• *r F /f 


(A.1 .5) 


where k Tf is the Thomas-Fermi screening wave vector and it is 
given by 




(A.1 .6) 


where f>(£) is the density of states including spin. 

The significance of k Tf r is not hard to find. If we have a 
bare Coulomb potential energy 

Vi?; - e Vr > (A - 1 - 7 ) 

then its Fourier transform is 

V (<p = Hwe 2 -/^- , (A.1.«) 


and the effective Fourier transform of the screened Coulomb 
potential is 

«... - 

l ^ * Tf= e 

By inversion, we can then show that this leads to the following 
effective potential energy 

z -fe TF r (A.1. 10) 

W(v) = j- e 

Many authors have shown for Bloch Wave functions 
referred to band n , that 


HTre t r l<£>i e ) 

6 ( “ *• V^ 2 - , C - fjO -t~t CO - l ex' 

4 4 M** M . c h,n-*- 


(A.1 .11) 
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For non parabolic but still spherical energy 
Fermi-Thomas wave vector is still given by Eq. (A.l 
and Dietl , 1978) with Ju 

(see Appendix A. 8). 


surfaces the 
.6) (Szymanska 

(A. 1.12) 
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Dielectric Function 


We have already discussed the dielectric function and given 
an expression for a semiconductor. Here we want to summarize the 
results for the band structure of Hg Cd Te. Similar material 
is discussed by Whitsell (1981) . 

It is usual to start by dividing the dielectric constant 
into two parts, one part from the ionic lattice via its 
polarizability (lat) and the other part from the electrons (el) 
due to inter and intra band effects. 


£i<£ju>) - £|ar ^ e| (A ' 2J ' 

The lattice part shows two mode behavior, but can be ana- 
lyzed with no particular difficulty and reference can be made to 
the literature. The electronic part is rather complicted but has 
been well summarized by Lehoczky, Broerman, Nelson and Whitsett 
(1974), Nelson, Broerman, Summers and Whitsett (1978) and 
Lehoczky, Szofran and Martin (1980). The later paper gives a 
summary of how screening and the Thomas-Fermi cruve vectors enter 
into the longitudinal-optical phonon scattering and in ionizing 
impurity and electron-hole scattering. Generally speaking we 




(A. 2. 2) 


where f b and f (q) , are chosen appropriately according ^ to the 
scattering mechanism (and notice with £ b n we have £ h C l = ) and, 

depending on the circumstances, more complicated frequency 
dependent forms may be necessary. 
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Appendix A. 3 
Friedel Sum Rule 


The Friedel sum rule relates phase shifts (see Joachain, 
1975 for a definition of phase shifts) of partial waves of 
angular momentum,? scattered from a spherically symmetric poten- 
tial to the charge of the scattering center. At large distances, 
the scattering center is screened and so the charge of this 
center also equals the magnitude of the charge of the total 
number of electrons which are attracted to or repelled from it. 

The Friedel sum rule as customarily written (Ziman, 1978) is 

■2- = */ir Ej ^r0^lb F ) , (A.3.D 

where kp is the Fermi wave vector and 2 is the charge of the 
scattering center (in units of the magnitude of the electron 
charge) . 

For the conduction electrons in a semiconductor we write 
(Stern, 1967) ; ^ 

(2-e-roJ ( <e -> xf dfr - (A - 3 - 2) 


where is the Fermi function. Assuming a sharp Fermi surface, 
the first equation follows from the second. 

The usual form for the screening length (Eq. A. 1.6) can be 
derived from the Friedel sum rule and the first Born approxima- 
tion. The Friedel sum rule is thought to be exact (Mahan 1981, 
p. 236). 
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Appendix A. 4 

The First and Second Born Approximation 


Here we are concerned with the Born series for the phase 
shifts of partial waves. The development of the appropriate 
mathematics can be obtained from standard references such as 
Joachain (1975) . The equations we need are 

= - kj U ( kY) Ydr, (A- 1 *- 1 ) 


and 


Z j(k/Y) pH,,) + r.r') 


(A. 4. 2) 




k ji (ku) '] x ikr > ) J 


(A. 4. 3) 


where r > and r < are the greater and lesser of r and r respec- 
tively and , rji are the spherical Bessel and Neuman functions. 

One iteration yields the first Born approximation for the 
phase shifts 


Sj° - - A 


V , 


(A. 4.4) 


(see Eq. 2) . 

Two iterations yield the second Born approximation 

(A. 4. 5) 


-La*, sr - -b~ 4 - 


(see Eq. 4) . 


Variational methods can yield even better approximations 
(Joachain, ch. 10, 1975). The following agrees with the Born 

series through second order. 




(A. 4. 6) 
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Appendix A. 5 

Band Stru cture of Ho, Cd Te 
The Three Lev el Hoael 


It is customary to neglect the bands not shown. 



The /f • Hamiltonian is solved with this model (Kane, 1957) 
near the T 7 point {k - 0) . The equation describing the conduction 
band (c) , the light hole (lh) band and the split off band is 


4- (b ~ * 


P 2 " k x ) £ — 3 hp fe. 


(A. 5.1) 


where p is a momentum matrix element. The heavy hole (hh) band 
is described by a simple parabolic band. _JThe spin-orbit interac- 
tion mixes mixed spin functions and the k * p interaction mixes s 
and p functions. The dispersion relations are non parabolic and 
the electron effective masses are typically small. 

Whitsett (1981) gives a two and a half page summary of 
Kane's band theory and discusses the band structure of the 
similar compound Hg, Cd Se. For X>X 0 = *' 7 (the cross over 
point) the band structure x is like InSb. and the light hole 
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part of V g are coupled by the interaction and this 

"inverts" tnem. For \ < X b , becomes a valence band and 
Tg(lk) becomes a conduction band. The degeneracy of these two 
below x o causes the energy gap to vanish. 
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Appendix A. 6 
The Frohlich Hamitonian 


In Hg 1 _ Cd Te, at not too low temperatures, the scattering 
of current 1 carriers by polar longitudinal optic phonons is 
particularly important. The Frohlich Hamitonian describes this 
interaction. As given by Mahan (1981, p. 487) it is 


H - - 2 w -< a - 


4- 


r 1 1 

z, P h’ c W<v ’ 


(A. 6.1) 


where 




are the customary creation and annilation 


operators for electrons and phonons and f(q) which is propor- 
tional to the effective charge for interaction (Ehrenreich, 1956) 
is defined by Mahan. This coupling is not effective for trans- 
verse modes because these do not create strong electric fields. 
It is usual to approximate to by two modes ( for HgTe and CdTe) . 

V 


In order to discuss screening properly, one needs to add to 
this Hamiltonian electron-electron interactions. Screening 
affects the way the lattice interacts and so it affects the 
phonon energies. Screening also affects the electron-phonon 
interactions (Enhrenreich, 1959) . Screening is included by 
insertion of the correct dielectric function (Zawadski, 1974, 
Lehoczky et al, 1974, Whisett, 1981 p. 98). 
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Appendix A. 7 
Approximations 


We first give an order of magnitude estimate for the 
quantity LaZ which we have assumed to be small. 


If L = 1A 


Z 

a 


Uni si (mkv 

fi 2 - 


m* = . 025in 

e 


then we find 


LaZ = .1 


Although L might be larger, we also should correct the above by 
dividing by a static dielectric constant so the estimate may not 
be too far off. If the electron has an energy of 1 electronvolt 
and a mass as given above we estimate (on a parabolic band 
assumption) that 

„ y u. -I 
£ ^ x / o ^ t 

thus if L is lX then 

= . 0? .. 

The following simple argument picks out what X s should be 
important. If b is the impact parameter we expect 

k L ~ -1. 

Now if b>L then there should be no scattering for those values of X 
for which 

b = */* >L, 

Thus if kL < I j 

we would only expect the £ = 0 terms to be important . 
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Joachain (1975 p. 172) notes that at zero energy 
Born series converges if 

f ^ 

Y lU(Y)lJr < I J ( X~°). 

O 

Substituting and evaluating the integrals we find this requires 

L 3 2- < /. 

Thus our basic approximations would appear to be O.K. 
provided kL and Lai are both small, which would at least be true 
at low temperature. 


( k = 0) the 

(A. 7.0 
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Appendix A. 8 
Density of States 
and 

Screening Length 


If f>l£) is the density of states per unit volume, 
authors give the screening length as 


L 


- z 
U) 



/ 


several 


(A. 8.1) 


where a background dielectric constant is included and 
Eq. A. 1.6 gives k TF (with A. 1.12). As we have seen, this can 
also be called the screening length in the 1st Born approxima- 
tion. This expression is still valid for non parabolic bands 
provided they are spherical. 


Now 




S . — ‘ijrlP'dk, 

(2 rr) A J 


(A. 8. 2) 


for spherical energy surfaces using the usual factor of (2*r) 
for the number of states per unit volume in real and ^ space and 
the factor of 2 comes from spin. This implies ' 






£ 4±l 

TT^ • 


(A. 8. 3) 


Sometimes the momentum effective mass (m*) is used. For m*. 


SO 


l 

U** ~ ^ dk 



if TT 2 -. 


(a. 8. A) 
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Appendix A. 9 
Better Approximations 


If we are to make better approximations we must eliminate 
the three approximations we made; namely, 

(a) m* in V(r) a constant for corrections to the 
first Born approximation 


(b) kL«l 

(c) only need ^ = 0 term. 

We first indicate what happens when we relax (a) and (c) . 
It is convenient to redo the derivation a bit to see how things 
fit together. 

We assume a m* is not a constant. In order to get 
agreement with previous results for the screening length, m* is 
interpreted as the momentum effective mass as in Eq. A. 8. 4. ir 
we assume . 

S, I =■ 

the Friedel sum rule can be written 

C * i 

s/tt 2 , (2trl) 


) Je = Z. (A.9.D 


With and B^, defined as previously 


S' 


4 





~t B ^ 



(A. 9. 2) 


If we let a = m*b where now m* is considered a 
putting the expression for ^ into the Friedel sum 


function of k, 
rule and using 


a' 

P 121 ' 


(A. 9. 3a) 


and 


(zi- 1)! '• 

Ax = ~ 


(A. 9. 3b) 
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(A. 9. 4) 


as long as \i\ « / , we find 



where as before L u > is the screening length in the first Born 
approximation. This is obviously not a simple equation to solve. 

To see what happens if we do not assume kL<<l, it is 
instructive to evaluate A 0 . 

flo = / jl Ut n i l Jr 


~ r/i_ , 

—r, (2a) re Jr. 

Uzy) 


(A. 9. 5) 


We easily obtain 


l - - | rfj A ( (+ 


(A. 9. 6) 


If kL«l we obtain, 

fi 0 =_ (Lai)l[/-3 (U) L 2 


(A. 9. 7) 
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which goes over to our previous results when (kL) is negligible 
compared to 1. Similar results can be obtained for , but they 
do not appear to be particularly transparent. a 
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